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ABSTRACT 

We report the discovery of 44.7 ms pulsations from the X-ray source CXOU J181335.1 — 174957using 
data obtained with the XMM-Newton Observatory. PSR J1813— 1749 lies near the center of the young 
radio supernova remnant G12.82— 0.02, which overlaps the compact TeV source HESS J1813— 178. 
This rotation-powered pulsar is the second most energetic in the Galaxy, with a spin-down luminosity 
of E = (6.8 ± 2.7) x 10 37 erg s . In the rotating dipole model, the surface dipole magnetic field 
strength is B s = (2.7 ± 0.6) x 10 12 G and the spin-down age r c = P/2P = 3.3 — 7.5 kyr, consistent 
with the location in the small, shell-type radio remnant. At an assumed distance of 4.7 kpc by 
association with an adjacent young stellar cluster, the efficiency of PSR J1813— 1749 in converting 
spin-down luminosity to radiation is » 0.03% for its 2—10 keV flux, w 0.1% for its 20 — 100 keV 
INTEGRAL flux, and « 0.07% for the > 200 GeV emission of HESS J1813-178, making it a likely 
power source for the latter. The nearby young stellar cluster is possibly the birthplace of the pulsar 
progenitor, as well as an additional source of seed photons for inverse Compton scattering to TeV 
energies. 

Subject headings: ISM: individual (HESS J1813-178, G12.82-0.02) pulsars: individual 
(PSR J1813— 1749) — supernova remnants 



1. INTRODUCTION 

The identification of high-energy (TeV) 7-ray emis- 
sion in our Galaxy with supernova products has opened 
up new prospects for understanding the energetics of 
these stellar remnants. Over half of the > 50 Galac- 
tic TeV sources are identified with supernova remnants 
(SNRs) or pulsar wind nebulae (PWNe) 1 . Several physi- 
cal mechanisms are available for generating the observed 
7-ray photons, including non-thermal bremsstrahlung, 
and inverse-Compton scattering of ambient microwave, 
IR, or optical photons off relativistic electrons. 7-rays 
can also be produced in hadronic collisions of high-energy 
protons off local material, from the decay of their neu- 
tral pion products. An unambiguous detection of the 
latter process associated with a Galactic SNR would pro- 
vide direct confirmation of cosmic-ray production in SNR 
shocks. 

A unique opportunity to explore these high en- 
ergy mechanisms is provided by the TeV source 
HESS J1813— 178, one of the brightest and most com- 
pact objects located by the HESS Galactic Plane Sur- 
vey (Aharonian et al. 2005a). Within the TeV ex- 
tent lies G12.82— 0.02, a previously uncatalogued young 
shell- type radio supernova remnant (Brogan et al. 2005; 
Hclfand et al. 2005). Furthermore, there is evidence for 
broad-band high-energy emission. Overlapping the SNR 
shell lies a 2-10 keV X-ray source AX J1813-178 (Bro- 
gan et al. 2005; Helfand et al. 2005), the INTEGRAL 
soft gamma-ray (20-100 keV) source IGR J18135-1751 
(Ubertini et al. 2005), and possibly GeV emission from 
the Fermi source 0FGL J1814.3-1739 (Abdo et al. 2009). 
Follow-up high resolution X-ray studies resolved the X- 

1 VHE 7-ray Sky Map and Source Catalog, 
http:/ /www.mppmu.mpg.de/" rwagner /sources/ 



ray emission into a point source and bright surrounding 
nebula, whose properties indicate that a young, ener- 
getic rotation-powered pulsar (Funk et al. 2007; Helfand 
et al. 2007) is responsible for the extended X-ray emission 
and, ultimately, the broad-band high-energy radiation. 
A deep radio timing search failed to detect pulsations 
(Helfand et al. 2007). 

In this Letter, we report on a long, continu- 
ous XMM-Newton X-ray timing observation of 
CXOU J181335. 1-174957, the central point source 
in G12.82— 0.02, resulting in the discovery of 
PSR J 18 13 -1749. We show that the spin-down 
luminosity of the pulsar is sufficient to power the 
broad-band X-ray and 7-ray emission. A nearby young 
stellar cluster (Messineo et al. 2009) may provide seed 
photons for upscattering to 7-rays, and is a possible 
birth place of the pulsar progenitor. 

2. XMM-Newton OBSERVATIONS AND RESULTS 

The Chandra source CXOU J181335. 1-174957 was ob- 
served by XMM-Newton with a 98 ks exposure on 2009 
March 27 using the European Photon Imaging Camera 
(EPIC; Turner et al. 2003). The EPIC pn CCD was 
operated in small-window mode (4'3 x 4f3 field-of-view; 
29% deadtime). This mode provides 5.7 ms time resolu- 
tion, allowing a search for even the most rapidly rotat- 
ing young pulsar. XMM-Newton is sensitive to X-rays in 
the nominal 0.1 — 12 keV range with energy resolution 
AE/E « 0.1/V-E(keV). The target was placed at the 
default EPIC pn focal plane location for a point source. 
Data was also acquired with the two EPIC MOS CCD 
cameras (MOS1 and MOS2). These were operated in full 
frame mode with a time resolution of 2.6 s, insufficient 
to search for a typical pulsar signal. The medium filter 
was used for both EPIC instruments. 
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TABLE 1 

Timing Parameters of PSR J1813-1749 

Parameter Value 

R.A. (J2000) a 18 h 13 m 35.166 s 

Decl. (J2000) a -17°49'57.48" 

Epoch (MJD) 54917 

Period, P (ms) 44.699297(4) 

Period derivative, P (1.5 ± 0.6) X 10" 13 

Characteristic age, t c (kyr) 3.3-7.5 

Spin-down luminosity, E (ergs" 1 ) . . . (6.8 ± 2.7) X 10 37 

Surface dipolc magnetic field, B s (G) (2.7 ± 0.6) X 10 12 

NOTE. — lcr uncertainties arc given. 
a Chandra ACIS-I position from Helfand et al. (2007). 
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Fig. 1.— Discovery of PSR J1813-1749 using XMM-Newton 
EPIC pn data of 2009 March 27. Contours of the signal power, 
realized by a Z\ search over frequency and frequency derivative 

(/,/), corresponding to the 68%, 90%, and 95% confidence levels 
for two interesting parameters. 



All data were processed using the SAS version xmm- 
sas_20060628_1801-7.0.0 pipeline, and were analyzed us- 
ing both the SAS and FTOOLS software packages. The 
observation was free of significant particle contamination 
and provided a near continuous 97.9 ks of good observ- 
ing time per detector. Photon arrival times from the 
pn CCD were transformed to the solar system barycen- 
tcr in Barycentric Dynamical Time (TDB) using the JPL 
DE200 cphemeris and the Chandra measured coordinates 
given in Helfand et al. (2007). 

To search for the expected pulsed signal in the pn CCD, 
2 — 10 keV photons were extracted from an 0.'5 diame- 
ter aperture centered on the source, optimized for the 
highly absorbed point-source embedded in the substan- 
tial nebular emission (see Helfand et al. 2007). A 2 25 -bin 
fast Fourier transform (FFT) search algorithm was used 
and a highly significant signal detected at P = 44.7 ms. 
We performed a refined Z\ (Rayleigh) search (Buccheri 
et al. 1983) centered on the FFT signal and localized 
the pulsed emission with peak power Z\ ~ 476. This 
corresponds to essentially zero false detection probabil- 
ity for the number of independent search trials. The 
pulse profile is slighly narrower then a sinusoid and shows 
no energy dependence in subdivided bands within the 
2—10 keV range. 

Given the long time span of the data, we also looked 
for the effect of a period derivative / on the pulsar signal 
by constructing a Z\ search over (/, /)-spacc. Figure 1 
shows the resulting contours for the signal power at the 
68%, 90%, and 95% confidence level for two interesting 
parameters. The formal result is / = 22.371716(2) Hz 
and ,/ = -(7.7 ± 3.1) x 10~ n Hz s _1 ; the uncertain- 
ties are lcr. This large / is detectable in the single 
long observation, as the quadratic term \f{t — t ) 2 in 
the phase ephemeris contributes —0.37 cycles of rota- 
tion over the elapsed time. However, its value is highly 



uncertain at the 95% confidence level, as shown in Fig- 
ure 1. A precise value of / will be determined eas- 
ily with a follow-up observation. The best value and 
lcr range for the spin-down parameters in the dipolc 
pulsar model are spin-down power E = —4ir 2 Iff = 
(6.8 ± 2.7) x 10 37 erg s _1 , surface dipole magnetic field 

strength B s = 3.2 x lO 19 ^-/// 3 = (2.7 ±0.6) x 10 12 G, 

and characteristic age r c = —f/2f = 3.3 — 7.5 kyr. The 
pulsar parameters are listed in Table 1. 

Figure 2 presents the pulse profile of PSR J1813-1749 
corresponding to the data folded at the peak signal 
value in (/, /)-space. The detected pulsed fraction 
for the extracted photons is 25 ± 3%, defined here as 
A(pulsed)/./V(total), where we choose the minimum of 
the folded light curves as the unpulsed level. The quoted 
uncertainties are derived by propagating the counting 
statistics of the light curve for 10 bins. 

The intrinsic pulsed fraction is difficult to deter- 
mine due to the background contamination from the 
PWN in the source aperture (see Figure 3). Using the 
high-resolution Chandra X-ray image of G12.82— 0.02 
(Helfand et al. 2007), we can resolve the pulsar emis- 
sion from that of the nebula for a given sized aper- 
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Fig. 2.— Pulse profile of PSR J1813-1749 folded in 20 phase 
bins at the frequency corresponding to the peak power in the (/, /) 
search pcriodogram shown in Figure 1. The dashed line indicates 
the combined contributions from the X-ray background, PWN, and 
detector contamination. Phase zero is arbitrary. 
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Fig. 3.— XMM-Newton EPIC MOS broad-band X-ray image 
of PSR J1813— 1749 and its wind nebula. The image is exposure 
corrected, smoothed, and scaled to highlight the PWN emission 
and faint sources; the pulsar (cross) is saturated in this display. 
Sources associated with the core of the young stellar cluster of 
Mcssinco ct al. (2009) are seen to the southwest (dashed circle). 



ture to estimate their relative contributions. Over the 
EPIC pn signal extraction aperture, allowing for the 
EPIC pn beam size, we derive a photon flux ratio of 
-Ppwn/.Fpsr+pwn = 0.51, after correcting for the Chan- 
dra background (« 15% of PWN counts). Assuming a 
similar relative spectral response between the two detec- 
tors, we use the Chandra ratio to correct the background- 
subtracted XMM EPIC pn counts to determine the PWN 
contamination in the source aperture. After correcting 
for the EPIC pn X-ray and detector background contri- 
butions, and taking into account the Chandra derived 
PWN contamination, we determine an intrinsic pulsed 
fraction for PSR J1813-1749 of 58 ± 9%. 

3. DISCUSSION 

3.1. Age and Energetics 

HESS J1813-178 is one of the more compact TeV 
sources to be associated with a SNR. The Gaussian 
extent of the source is only a — 2.2' ± 0.4' (Aharo- 
nian et al. 2006a), while the radio shell of G12.82-0.02 
is only 2.5' in diameter (Brogan ct al. 2005). It is 
not clear if these are significantly different. Unlike 
PSR J1813-1749/HESS J1813-178, many TeV sources 
associated with pulsars show a notable displacement be- 
tween the two, e.g., PSR B0833-45/Vela X (Aharo- 
nian et al. 2006b), PSR B1823-13/HESS J1825-137 
(Aharonian et al. 2006c; Pavlov et al. 2008), PSR 
J1809-1917/HESS J1809-193 (Aharonian et al. 2007; 
Kargaltsev & Pavlov 2007), PSR J1718-3825/HESS 
J1718-385 (Aharonian et al. 2007; Hinton et al. 2007), 
and possibly PSR J1617-5055/HESS J1616-508 (Aha- 
ronian et al. 2006a; Landi et al. 2007). For these sys- 
tems, the TeV source is larger than the X-ray PWN. 
Apparently, the electrons responsible for TeV emission 
via inverse- Compton scattering have longer lifetimes, and 
may have been injected at earlier times, than the X-ray 



synchrotron emitting electrons now in the inner PWN. 
It is also possible that the PWN is displaced in one di- 
rection by an asymmetric reverse shock resulting from a 
supernova that exploded initially in an inhomogeneous 
medium (Blondin et al. 2001). These systems contain 
"middle-aged" pulsars, > 10 4 yr, which are old enough 
for the reverse shock to have had its effect, or for rela- 
tivistic electrons to have diffused out of the high B-ficld 
PWN. Theoretically, a TeV source powered by a pulsar 
can live for up to 100 kyr as an inverse Compton emit- 
ter in regions of low (/iG) -B-fields, with luminosity even 
exceeding the current spin-down power of the pulsar (de 
Jager & Djannati-Ata'i 2008). Their TeV luminosities 
can then exceed their X-ray luminosities by a large fac- 
tor, which is observed for several sources. 

In the case of PSR J1813-1749/HESS J1813-178, the 
compact TeV source and lack of relative offset can now be 
explained by its youth, like the Crab and several young 
pulsars with compact HESS sources that are co-located 
with their X-ray PWNe. The other compact sources are 
PSR J1846-0258/HESS J1846-029 in the SNR Kes 75 
(Djannati-Atai et al. 2008), PSR J1833-1036/HESS 
J1833-105 in G21. 5-0.9 (Djannati-Atai et al. 2008), 
and the newly discovered 52 ms PSR J1747-2809/HESS 
J1747-281 in GO.9+0.1 (Camilo et al. 2009; Aharonian 
et al. 2005b). PSR J1846-0258 is a 0.325 s pulsar with 
a characteristic age r c of only 728 yr, and a spin-down 
luminosity E = 8 x 10 36 erg s" 1 (Gotthclf ct al. 2000). 
PSR J1833-1036 has r c = 4900 yr and E = 3.3 x 10 37 
erg s -1 . However, the actual age of PSR J1833— 1036 
is probably only ~ 870 yr as measured by the expan- 
sion of G21. 5-0.9 (Bietenholz & Bartel 2008). PSR 
J1747-2809 has E = 4.3 x 10 37 erg s" 1 and r c = 5300 yr. 
PSR J1813— 1749 shares with these young pulsars a char- 
acteristic age ^ 5000 yr, but its true age may be signif- 
icantly smaller than this because its initial period Pq is 
not necessarily << P. In fact, Brogan ct al. (2005) es- 
timate an age of < 2500 yr for G12.82-0.02 based on 
its size, and possibly as small as 300 yr if it is still in 
the free-expansion stage. In these relatively young and 
compact PWNe the high -B-fields make for efficient syn- 
chrotron X-ray sources, and inefficient inverse-Compton 
TeV emission. The Crab, PSR J1833-1036, and PSR 
J1846-0258, have ratios L x /L 1 of « 120, 40, and 10, re- 
spectively (Djannati-Atai et al. 2008). This ratio is « 0.5 
for PSR J1813-1749. 

The efficiency of PSR J1813— 1749 converting spin- 
down luminosity to radiation at a distance of d = 4.7 kpc 
is « 0.03% for the 2-10 keV X-ray flux, w 0.1% for the 
20-100 keV INTEGRAL flux, and w 0.07% for the > 200 
GeV emission of HESS J1813-178 (see Table 2). Mat- 
tana et al. (2009) used empirical correlations of the ratio 
L x /L 7 with E among known HESS PWNe to estimate 

that E ~ 1.5 x 10 37 erg s _1 and t c <~ 6 kyr for the as-yet 
undiscovered PSR J1813— 1749, which turns out to be a 
reasonably accurate prediction of their actual values. 

3.2. Environment and Associations 

The location of G12.82— 0.02 near a young stellar clus- 
ter rich in massive binaries and containing a second su- 
pernova remnant (G12.72-0.00) suggests an association 
(Messineo et al. 2009). Those authors determined clus- 
ter membership for giants and early-type stars, and de- 
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TABLE 2 

Multi-waveband Luminosity Measurements of PSR J1813-1749 



Mission 


Energy Band 


Source 


Luminosity 3, 


£ b 


Reference 








(erg s" 1 ) 


(%) 




Chandra 


2-10 keV 


AX J1813-178 


1.9 x 10 34 


0.03 


Helfand et al. (2007) 


INTEGRAL 


20-100 kcV 


IGR J18135-1751 


7.3 x 10 34 


0.1 


Dean & Hill (2008) 


Fermi c 


0.1-1 GeV 


OFGL J1814.3-1739 


5.4 x 10 35 


0.8 


Abdo ct al. (2009) 


HESS 


> 200 GcV 


HESS J1813-178 


4.8 x 10 34 


0.07 


Aharonian et al. (2006a) 



a Total luminosity in the band assuming d = 4.7 kpc. 

b Emciency in the band assuming spin-down power E = 6.8 X 10 37 erg s — 1 . 

c Conversion from Fermi photon flux assumes a power-law photon index T = 2. How- 
ever, association of the Fermi source with G12.82-0.02/PSR J1813-1749 is not yet 
established (see §3.2). 



rived a consistent distance of 4.7 kpc using radial ve- 
locities and optical/infrared extinction. They estimated 
an age of 6-8 Myr, and a total initial cluster mass of 
2000 — 65OOM0. It is plausible that the progenitors 
of both G12. 82-0.02 and G12. 72-0.00 were born in the 
cluster, and had masses similar to those of the red su- 
pergiant and WR stars now present, 20 — 30M Q . The 
cluster is centered 4.4' southwest of G12.82— 0.02 (see 
Figure 3). At this offset, a 5 kyr old neutron star born 
in the center of the cluster requires a transverse velocity 
of w 1200 km s _1 to reach its present location, or less if 
it was born in the outskirts close to its present location. 

It is plausible that seed IR/optical photons 
from this star cluster enhance TeV emission from 
HESS J1813— 178, via inverse Compton scattering in 
its extended PWN, above the minimum expected from 
the cosmic microwave background. Such a role was also 
hypothesized by Helfand et al. (2007) for the adjacent 
W33 star-forming region, at an estimated distance of 
4.3 kpc. Although these objects are not necessarily all 
at the same distance, they could be. 

Pulsed emission from PSR J1813— 1749 is expected 
to be detectable by Fermi, and the GeV band 
could be where the luminosity peaks, as it does in 
many energetic pulsars. Although a possible asso- 
ciation of the Fermi source OFGL J1814. 3-1739 with 
HESS J1813-178/ G12.82-0.02 has been noted by Abdo 
et al. (2009), the SNR lies just outside the large, 11' ra- 
dius 95% confidence Fermi error circle. It is not yet 
possible to identify definitively the GeV emission with 
the pulsar, which would require detection of pulsed 7- 
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rays. The implied luminosity of OFGL J1814.3-1739 (as- 
sumed isotropic) in Table 2 is consistent with the trend 
of efficiencies among 7-ray pulsars (Thompson 2004). 
Accumulation of more exposure with Fermi and detailed 
investigation of this region could be informative. 

4. CONCLUSIONS 

Using XMM-Newton, we discovered PSR J1813-1749, 
a highly energetic 44.7 ms pulsar associated with the su- 
pernova remnant G12.82— 0.02 and powering its PWN 
from X-ray to TeV energies. Its preliminary spin-down 
properties make it one of the most energetic pulsars in 
the Galaxy, possibly second only to the Crab. Its appar- 
ent youth explains the small extent of its SNR and TeV 
nebula relative to older pulsars. The high E and young 
age support a one-zone synchrotron/IC model of the type 
applied by Funk et al. (2007), in which PSR J1813-1749 
accelerates electrons to > 10 15 eV. The lack of an X- 
ray SNR shell, coupled with a high spin-down power for 
PSR J1813-1749, argues that the PWN, not the SNR, 
is almost certainly the TeV source. 
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